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Predicted by Wannier in 1960 band states quantization in a constant electric field, E n = const 
±n£", where n — 0, 1, 2, and £ is proportional to the strength of electric field [this kind of 
spectrum is commonly referred as the Wannier-Stark ladder (WSL)], implies that the probability 
of tunneling through a tilted band should have £ spaced peaks, at least, under the weak coupling 
of the band states to the source and drain electrodes. It has been shown, however, (Phys. Rev. 
B 63, 2001), that the appearance of the canonical WSL is preceded by WSLs with other level 
spacing, namely, £ m i / m /(l — 2m'/m), where m and m' < m/2 are positive integers specifying certain 
applied voltage. Here we show that canonical and noncanonical WSLs, in addition to different peak 
spacing in the transmission spectrum, have other pronouncedly distinctive features. As an example, 
for the former, the peak and valley tunneling probability decays exponentially with the increase of 
applied voltage. The corresponding exponents are given by the sum and difference of two Fowler- 
Nordheim-type exponents implying an anomalous increase of the peak-to-valley ratio. These and 
other signatures of extended states (es), surface localized states (sis), and Wannier-Stark (WS) states 
in the through tilted band transmission spectrum are discussed on the basis of (derived for the first 
time) nearly exact explicit expressions of es-, sis-, and WS states assisted tunneling probability. 
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I. INTRODUCTION 



Introduced by Wannier in 1960,H'EI the concept of Wannier-Stark ladder (WSL), associated with the spectrum 
of electron band states in a constant electric field, gave rise to an unobservable literature of discussions aimed at 
the development of in-depth understanding of electric field effects in solids. In nineties, this activity has resulted in 
several repopjts,on experimental observation of the Wannier-Stark (WS) effect, Bloch oscillations, and related coherent 
phenomena,0~Ei appearance-of new ideas (such as, e.g., dynamic localization^ and quasi-energy band collapse in time- 
dependeat .electric fields Jxla) and new possibilities to study static and time dependent field effects in optically driven 
lattices-tsHta 

In its canonical form, the WSL represents a spectrum of evenly spaced electronic levels E n = const ±n£ , where 
n = 0, 1, 2, and £ (the level spacing proportional to the strength of electric field) is equal to the Plank constant 
times Bloch oscillation frequency. In the tight-binding description, this kind of spectrum is readily expectedJiS if the 
constant shift £ of electron on-site (atomic) energy along the electric field is comparable with the width of the parent 
zero-field band of Bloch electron states. Such an extreme case of WSL, though is easy to think about, is of little, if 
any, physical interest, since in most cases, the band state spectrum emerges as a result of delicate interplay of electric 
field and interatomic interaction effects, where the methodology of perturbation theory is inappropriate. 

In early seventies, it was realized that the tight-binding model presents a unique opportunity to treat field effects 
in solidsj-on—rigorous basis, since its Hamiltonian admits obtaining the exact formal solution of the Schrodinger 
equatior£j~E9 and hence, the Green function prpblera f| rp without any restriction on the field parameter. A number 
of useful results has been obtained-ui this wayJi3iia23iHil in particular, that in finite systems, some minimal voltage 
is required to open the WS band.Ej Unlike the Wannier model, which refers to the bulk electron states only, the 
field-induced surface states also have been brcaight. tp light. Itjaps shown that eigenenergies near the spectrum edges 
behave similarly to zeros of the Airy functionEalflo Recently,c3 within the same model, new regularities for the bulk 
and surface electron states influenced by a constant electric field were predicted. The bulk levels, for instance, have 
been shown to form (under certain voltages) WSLs with noncanonical level spacing q£, where q (> 1) can be an 
integer, as well as a fractional number. Thus the concept of WSL appeared even more rich and fertile physically, than 
was commonly thought. 

Next step to be undertaken is to specify such physical quantities, which allow one to observe distinctions between 
the true Wannier levels and the others, for instance, the levels that correspond to extended states (es) and form 
noncanonical WSLs, and those corresponding to surface localized states (sis), whose spacing is governed by three 
different lawsBa Here, in focus is the probability of tunneling (or transmission) of electrons through a tilted tight- 
binding band. This quantity is of relevance to the phenomenon of electrical breakdown of ultra-thin dielectric layers 
and a number of other electric field effects in solids. The energy dependence of electron transmission can be directly 
probed by means of the ballistic-electron-emission microscopy technique.Ej And properties of through tilted band 
tunneling are intimately related to Zener tunneling and Franz-Kcldysh effect. 

II. MODEL AND GENERAL CONSIDERATION 

We choose to discuss the electric field effects on tunneling through a spatially finite and energetically restricted band 
of electronic states in the context of electron ballistic transport in superlattices. The profile of TV-well supperlattice 
model potential is shown by dashed line in Fig. 1. Due to electron confinement in a well, the well spectrum is 
quantized. Due to tunneling between neighboring wells, the quasi-discrete (in its lower part) spectrum of an isolated 
well splits into minibands of allowed electron energies, separated by the energy gaps, which do not contain electron 
states. Since the hight and width of barriers can be engineered freely within a wide range, it is well possible to have 
the lowest miniband separated from the others by an energy interval much exceeding its band width E^. An example 
of this kind is represented in Fig. 1. 

Under a restriction that the electrostatic potential drop across the superlattice does not exceed much the magnitude 
of £\,w) the description of electron transport along the superlattice can be started from the Hamiltonian matrix for a 
single band 

H nn , = £ [n - (/V + l)/2] S n , n > + *|n-n'|,l, (!) 

where indices n, n' — l,J\f number the wells in the superlattice (see Fig. 1); £ — eFa/ f3; e, F, and a are, respectively, 
the absolute value of electron charge, electric field strength, and periodicity of the superlattice a = w w + Wb! the 
energy of electron resonance transfer between neighboring wells /3 serves as an energy unit; the zero- field energy of 
the lowest level in the non-interacting wells (the electron on-site energy) is set equal to zero. 
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The eigenvalues of the Hamiltonian matrix (1) are lying within the (dimensionless) energy interval (along E axis 
in Fig. 1) which is equal to E^ w + eV, E^ w — 4, eV = £{N ' — 1). Excluding classically inaccessible regions in the n-E 
'coordinates', one obtains a tilted band of electron states. Such a band is shown in Fig. 1 for two voltages eV < E hw 
and eV > E^. Sloped lines on the band diagrams indicate the n-E (unshaded) regions, where the probability to find 
an electron with the given energy in a well n (= 1, 2, Af) is proportional to the tailing part of the corresponding 
squared eigenfunction of H nn i . As is seen from the figure, if the applied potential is smaller, than the miniband 
width in zero field, eV < £W, the band spectrum can be divided into es-band and two sis-bands; at higher voltages 
eV > Ehvr, the es-band is replaced by WS-band. 

Let us assume that from its opposite boundary layers, the superlattice is contacted homogeneously to ideal semi- 
infinite drain and source electrodes and, that the interaction between its first (7Vth) well and the drain (source) lead 
can be described by a single parameter V. According to the model, the electrostatic potential energy inside the leads 
is constant (though different) and drops from the source to drain by value of eV , entirely in the region occupied 
by the superlattice, see energy diagrams in Fig. 2. Suppose also that the spatial variables of the Hamiltonian of 
the system drain-superlattice-source are separable. Then, the electron energy E can be broken into a sum of two 
non-mixing components which correspond to the energy of electron motion parallel and perpendicular (E±) to 
the electric field, the latter directed along the superlattice. In such a system, an electron wave, having the energy E 
and propagating freely, say, in the source towards the draicij-spll be transmitted through the superlattice with the 



probability T{E,£). Due to the Landauer-Biittiker thearyJiJcj the transmission probability is directly related to the 
current- voltage relation, see, e.g., the Datta overview.EJ 

Probably, for the first time, the model of ballistic electron transport outlined above was suggested and solved on the 
basis of the Green function technique by Caroli et a/.Ell Originally, it was applied to the description of tunnel current 
in metal-insulator-metal heterostructures. To calculate the current, the Green function has to be found. In the cited 
work, however, the latter was not specified. Later on, the variations of Caroli et al. treatment have reemerged in 
a nunabeiijof physical contexts. In particular, it has been used to examine the quantum conductance of molecular 
wires. ErH 

For the model at hand, the transmission probability T{E 1 £) is conveniently expressible in terms of certain matrix 
elements of the Green functions referred to the non-interacting source and drain leads, and the scattering region, 
which here is the superlattice. The miniband width is normally much smaller, than that of typical contacting leads. 
Therefore, in the actual energy interval, the variation of the lead-related Green functions can be regarded as negligible 
in most cases of interest. Thus, the only energy dependent quantities remaining in the transmission probability are 
those referred to the superlattice, i.e., to the Hamiltonian (1). 

Combined together, the model assumptions yieldEJ 



4A 2 Gj^(E) 

[1 + A 2 Ga(E)] 2 + A^G^E) - G^ME)} 2 + ^G\ N {E) ' 



where Ga(E) — G\ : i(E)Gjsf : j^(E) — G\ N -(E); and the (superlattice) Green function obeys the equation 

N 

^ (E8 ni n» — H nn ") G n " n '(E) — 5 ntn >. (3) 
n"=l 

In Eq. (2), the quantity A ~ V 2 (independent of energy) associates with the weighted local density of states on 
the drain/source-superlattice interfacial surfaces. The Green function matrix elements, which appear in the above 
equations, depend on E\\ and £ . These are the only energy dependent quantities in Eq. (2) and therefore, here and 
henceforth the subscript in Eu , as well as an indication of the Green function xLependence on £, are omitted for brevity. 
Analytical expressions of Gn(E) = —G^fj^(—E), and G\m(E) are given byEil 

Djs(E,£)G n (E) = J t i+(AT-i)/2(z)Y fl -(M+i)/2(z) - 5^+(jV-i)/2(«)J A t-(AT+i)/2(«) ) ( 4 ) 



Dat{E 7 £) = J^+(Ar+i)/2(z)Y fi _ {A f +1)/2 (z) - Y^ +{Ar+1)/2 (z)J fl _ {A f +1)/2 (z), (5) 

and Dj^{E,£)G\jsf(E) = £/ir, where /i = E/£, z = 2/£, and J M (z) and Y fl (z) are the Bessel functions of the first and 
second kind, respectively. 

Before going into a detailed discussion of T(E 7 £) dependence on energy, field strength, and superlattice length 
(the latter is not indicated explicitly), some relevant remarks might be useful. As is mentioned above, the n-E 
(length-energy) areas, which are classically accessible for electrons (shaded in Figs. 1 and 2) are distinct, depending 
on whether the electrostatic energy is smaller or larger, than the zero-field band width. For this reason, the voltages 
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restricted by the conditions eV < E^ w and eV > E\, w will be referred as low and high voltages, respectively. The 
qualitative difference between the two cases stems from that the electron states extended over the entire superlattice 
(indicated in the figures by 'es') can exist at low voltages, but not at high voltages, when the es-band is shrunken and 
the WS-band opens instead. 

Contrary to extended states, which "connect" the source and drain electrodes, the electron states located within 
either of two triangular areas, are separated from the other side_of the superlattice by a triangular barrier, see Fig. 2. 
(The field induced surface localized states exist if £ > 12.2/A/" 3 .l3) Two triangular barriers separate the states located 
within the WS-band. Hence, depending on the applied voltage and energy of electrons to be transmitted through the 
superlattice, the electron flux can encounter no barriers at all, one triangular barrier, two triangular barriers, or a 
trapezoid barrier, as is exemplified by arrows labeled in Fig. 2 by RT, FN, TB, and WS. The through trapezoid barrier 
(TB) tunneling, which-also includes the Fowler-Nordheim (FN) regime of tunneling, has been extensively studied in the 
WKB approximation.^ Such an approach is inapplicable for the description of the present discrete, finite band model. 
In particular, it fails to reproduce the resonance structure of T(E,£), which might be expected for the sis-assisted 
tunneling (i.e., tunneling through a triangular well). For the resonance tunnelins-(RT) transmission, the energy and 
length dependence of T(E,£) in zero field £ = has been analyzed previously!^! However, the presence of electric 
field gives rise to new regularities which have not been studied so far. In what follows, the transmission spectrum will 
be examined in the RT, FN, WS, and TB energy intervals, where surprisingly accurate explicit expressions of T(E,£) 
can be derived from the model exact equations (2), (4), and (5). 

III. EXTENDED STATES ASSISTED TUNNELING 

At low voltages the mid band levels can form WSLs with noncanonical level spacin; 

El = {l + l/z} l-2m>/m £m '/ m > (6) 

where the upper (lower) factor should be used for an odd M = 2N + 1 (even M = 2N) number of wells in the 
superlattice; I = 0, 1, 2, m and m! < m/2 are positive integers, which specify the value of the field parameter 

(Af-l)£ roVro = 4cos(7r— ) . (7) 

For some voltages, the correspondence between the approximate WSL energies and exact eigenvalues (solutions to 
T>j\f(E,£) = 0) can be seen in Fig. 3, where the former and the latter are indicated by circles and stars, respectively. 
In our example, Af — 101; the increase/decrease of TVlimproves/worsens the accuracy of Eq. (6). 

As is implied by the derivation of Eqs. (6) and (7)j£j the WSL with the level spacing, equal to £ times any but > 1 
integer or ratio of two integers, appears in the spectrum mid under a voltage, determined by Eq. (7), within the range 
< eV < £"bw This result may be thought as a reminiscence of approximately equidistant levels with the spacing 
equal to ir/N near the center of zero-field band. For an experimental observation, the WSLs of the type Ei = lq£, 
where q is equal to either a positive integer > 2 or 1 + 2m'/ (m — 2m') with a not too small second summand, are 
perhaps most meaningful, because the level spacing is well distinguishable from that in the canonical WSL. It also 
may be of importance from the experimental point of view that any of (noncanonical) WSLs can be obtained only at 
a unique voltage. For instance, the WSL with triple-f level spacing should appear, when the electrostatic energy is 
exactly equal to E^/2, if N is odd, or to E^ w /2 — £, if N is even. As we believe, these results can help to identify 
the electric field effects. 

To obtain an explicit expression of T(E,£), we use in Eq. (2) an approximation of the Green function matrix 
elements at the energies E n = n£ with n equal to a positive integer restricted by tha condition n < < Af. Performing 
some algebra, which is much in the spirit of the derivation of noncanonical WSLsjH3 we get instead of Eqs. (4) and 
(5) 

Gn(£/ n ,£)« r-77. — v — > (oa) 

sm(2nx) 

G^(E n ,£)^- S J^±M, (8b) 
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(8c) 



where \ — arccos(£iV/2). Unlike the exact expressions for the Green function matrix elements, Eqs. (8a), (8b), and 
(8c) make sense only for large values of M (set to be odd), and only for the indicated energies. Using the above 
expressions in Eq. (2) we obtain 



4A 2 sin 2 x 

A A 2 sin 2 x + sin 2 (2nx) [1 + 2 cos(2x)A 2 + A 4 ] ' 



T(E = n£,£)^ 77¥ - r - 2 , 2 . o __ AM n _^, A9 , HT > (9a) 



or, equivalently, 



T(n£ £) „ A 2 (A-£ 2 N 2 ) 

(U ' J ~ A 2 (A -£ 2 N 2 )+ sin 2 [2na.rccoa(£N/2)}[(l~ A 2 ) 2 + A 2 £ 2 N 2 }' ( ' 



Note that Eq. (9a) is also valid for high voltages, in which case x = 

cosh" 1 (£'A^/2), and sinx and cosx should be 

replaced by sinhx and coshx, respectively. 

According to Eq. (9b), if the applied potential satisfies Eq. (7), the transmission coefficient is equal to unit for 
E = n£ with n divisible by m, if m is odd, or by m/2, if it is even. This result makes it obvious that T{E 1 £) has a 
resonance-like structure, as a function of the applied voltage at a fixed energy and, as a function of energy at a fixed 
voltage. Thereby, noncanonical WSLs are exposed in the transmission spectrum, see Fig. 3, as almost equidistant 
peaks. It is seen that the level spacing is in a good agreement with formula (6) (the WSL energies are indicated in 
Fig. 3 by circles on the E-axis). However, the above approximation does not reproduce the decrease of transmission 
peak maxima, which is a pronounced tendency of the exact dependence of the transmission probability on energy. 
Represented by a few graphs in Fig. 3, it reflects a general trend of resonance tunneling phenomena, which implies the 
unit transmission probability only for totally symmetric systems (e.g., as conventional barrier- well-barrier symmetric 
heterostructures) . In the given case, the system symmetry is broken by the applied potential, the increase of which 
results in an increasing suppression of the transmission peaks close to the es-band edges. 

With the increase of superlattice-to-lead coupling parameter \ A\, transmission peaks broaden, see Fig. 3, which is 
a kind of expected behavior from other resonance tunneling structures. Figure 3 also exposes intuitively not expected 
behavior of the transmission spectrum response to increasing voltage. As is clearly seen from the comparison of the 
peak positions for ever lager potential (from bottom to top), the peak spacing decreases with the increase of £. So, in 
addition to a distinctive level spacing, and hence, peak spacing, which is characteristic for noncanonical WSLs, they 
are also distinguishable, from the experimental point of view, due to a specific field dependence of the transmission 
peak spacing. The latter decreases with eV, while in a canonical WSL it should increase. 



IV. SURFACE LOCALIZED STATES ASSISTED TUNNELING 

By using the standard approximations of Bessel functions with large arguments and small or large orders it can be 
shown that the exact expressions of Gn(E), Gj^j^(E), and G\j^{E) [see Eqs. (4) and (5)] are accurately reproduced 
within the sis-band energy interval |2 — eV/2\ + £ < E < 2 + eU/2 — £ by the following relations 



V N {E,£)G X1 (E,£) 
V M {E,£)G NM (E,£) 



£/ir 



\/sin £ sinh <5 



exp ( -$ 5 ) { 



f 2 ^ * 

cos {-* e - + 

cos {-* t - 

'2 , 7T 



(10) 



where 2 cosh 5 = E + eV/2, $a = S cosh 6 - sinh 6 , 2cos£ = E — eV/2 (0 < £ < n), and $ 5 = sin£-£cos£. Note that 
the excluded £ intervals can contain not more, than one sis level each. 

Using the above expressions for the Green function matrix elements in Eq. (2), we get 



T(E,£) 



A A 2 sinh 6 sin £ exp {—26) 



(1 + A 2 e- 2S ) 



7T 

4 



A 2 < 



exp 



(11) 
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showing a resonance-like structure modulated by a function exp (— 4$g/£) that decays exponentially with the increase 
of energy. 

For the energies of sis-band levels Ef s , eiven by solutions to equation T>_\f(E, £) ~ cos (2$^/£ — ir/4 + £) =0 within 
the energy interval |2 — eV/2\, 2 + eV/2fB expression (11) simplifies to 

-.sis c \ _. 4sinh(5 n exp(-2<5„) _( 4_ 



sm£„ [1 + A 2 exp{-26 n )\ \ £ J 

where S n = cosh -1 + eV/4) and = arccos (£^ ls /2 — eV/4) . In the case of weak superlattice-to-lead coupling 

A 2 « 1, Eq. (12) determines local maxima of T(E,£). Thus with the replacement E^ s — > E, the r.h.s. of Eq. (12) 
gives a function enveloping peaks of the transmission spectrum in the region of sis-assisted tunneling. 

As is seen in Fig. 4, except a small £ interval above the top of es- or WS-band, an explicit analytical description 
of the sis-assisted tunneling given by Eq. (11) is indistinguishable from exact calculations. This is a central result of 
the section. It might be worth emphasizing that details of the resonance structure of T(E,£), such as the position 
of peaks, and their width and intensity, are model dependent. In contrast, factor exp(— A$>$/£), which prescribes for 
the sis-assisted tunneling probability an exponential decrease with the increase of energy (counted from the spectrum 
center) and an exponential increase with the increase of electric field strength, is characteristic for the linear drop of 
the electrostatic potential within the scattering region and does not depend on a particular model of the interface 
(connection) with the leads. The reason for this is that sis-assisted tunneling is controlled by exponential tailing of 
the sis wave function in the classically forbidden region. Therefore, the exponential decay of the through sis-band 
tunneling probability can be evaluated from the ratio of the probabilities to find an electron with, say, the energy E^ s 
in the first and last wells of the isolated supperlattice. Such an approach, which is much more simple technically (since 
the leads and connection to thern—are out of consideration), gives a reasonable result even for the pre-exponential 
factor of the enveloping functionO T(E^ S ,£) = (sinh<5„/ sin£ ra ) exp(— 2S n ) exp (— 4$a ri /£), which differs from the 
correct expression (12) only by factor of four. 

To make Eq. (11) easy readable, one can use an approximate expression 

3 ® S ~\(E + eV/2)y 2 , eV>E bwi (13) 

where eV = \eV — E bw \ = \eV — 4| is the excess of electrostatic potential energy over zero-field band width, and eV/2 
has the meaning of the top of es-band (WS-band) for low (high) voltages. The above approximation reproduces the 
exact dependence reasonably well up to eV of the order of -Ebw, see Fig. 5. 

Using Eq. (13) in Eq. (11), the exponential factor exp (— 4$^/£) can be replaced, if eV < E bw , by. 
exp [— 4(E — eV /2) 3 / 2 /(3£)] , which is identical to that appears in the Fowler-Nordheim theory of field emission.L3 
At high voltages, however, factor exp [~4{E + eV) 3 / 2 /{3£)] does not have a semi-classical analogy. In some more 
details, links of Eq. (11) with the probability of tunneling through a triangular barrier are discussed in Ref. 37. 



A. Effects of coupling on sis transmission spectrum 



The superlattice-to-lead connection (represented in Eq. (11) by a single parameter A) is dependent on a number 
of factors. For this reason, in relevant heterostructures the magnitude of effective coupling \A\ may vary by orders. 
It is of interest therefore, to trace the dependence of the sis-band transmission spectrum on the coupling strength. 

A common expectation is that with the increase of coupling strength the resonance structure is smeared out. This 
is really true for resonance tunneling in zero field. By contrast, as can be readily seen from Eq. (11), such an 
expectation is not justified for sis-assisted tunneling. In the latter case, the peak position and their sharpness is 
essentially determined by two cosine terms in the denominator of r.h.s. in Eq. (11). In the absence of anyone of the 
two, the transmission spectrum would contain infinitely high resonances either at the sis energies E^ s (A = 0, the 
extreme case of weak coupling) or at zeros of cos (2$^/£ — tt/4) (the extreme case of strong coupling). 

If the coupling is weak, the peak spacing repeats the level spacing within sis-band. In due course, the latter is 
proportionaL-tp £, 2£, and ruled by poles of the Airy function, within the lower, mid, and upper part of the sis-band, 
respectively.^ These regularities can thus be observed in the sis-band transmission spectrum, as is illustrated in Fig. 
4a by open circles {£ spacing), filled circles (2£ spacing), and squares (Airy type spectrum). 

In the case of strong coupling A 2 >> 1, the position of peaks in T(E,£) is determined by zeros of the second 
summand of the denominator in Eq. (11), and we have instead of Eq. (12) 
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4sinh<L / 4 \ . . 

T(E P ,S) = — ^ eX p(--^j, (14) 

where the values of £J P are given by solutions to equation cos (2$^/£ — 7r/4) = 0. These are shifted with respect 
to sis energies £^ ls , and can be shown to obey the same regularities as those, observed for the eigenvalues. The 
hight of peaks does not differ much in the cases of strong and weak coupling, as can be seen from comparison of 
Eqs. (12) and (14). At the same time, if \A\ is large, the wells between peaks are approximately A 2 times deeper. 
Summarizing, in the case of sis-assisted tunneling, strong coupling with the leads makes the resonance structure of the 
transmission spectrum even more pronounced. Otherwise, the weak and strong coupling result in a similar structure 
of the transmission spectrum with identical regularities of the peak spacing dependence on the electric field, sec Fig. 
6. 

On the other hand, it follows from Eq. (11) that the resonance structure will be essentially smeared out in the 
case of intermediated coupling strength, when \A\ sa 1. This unusual behavior of sis-band transmission spectrum is 
illustrated by calculations of T(E,£) for different couplings, \A\ « 1, \A\ = 1, and \A\ >> 1 in Fig. 4b. One can 
see that with the increase of coupling strength, the transmission spectrum, at first, losses its resonance structure, and 
then acquires it again with, roughly, interplaced peaks and wells, and deepened wells. 

The above analysis shows that the resonance tunneling assisted by surface localized states, which appear as a result 
of the band tilting by the applied voltage, is characterized by a kind of unique (at least, met not often) dependence 
of the transmission spectrum on the coupling with electron reservoirs from and to which the tunneling occurs. 

Some minimal voltage is required for the first sis to appear and hence, sis-assisted tunneling to be possible. Further 
increase of the applied potential results in sis-band opening up to its maximal width i?bw ■ The increase of sis-band 
width is accompanied by the appearance of peaks in T(E, £) with the spacing governed by the following regularities. 
For eV < Ebw/2, the peak spacing is close to that of Airy spectrum. For E^/2 < eV < £bw, the Airy type 
peak spacing gradually changes to the double £ peak spacing characteristic for H nn / eigenvalues in the mid of the 
maximal- width sis-band. Finally, for eV > E^ Wl there exists the third characteristic energy interval (close to the top 
of WS-band), where peaks of the transmission probability are £-spaced. The latter spacing is commonly regarded as 
the WSL trademark. 



V. WANNIER-STARK STATES ASSISTED TUNNELING 



In the energy interval of bulk states (i.e., in the mid of the full spectrum), switching from the low to high voltages 
results even in a more profound restructuring of the transmission spectrum. This can be expected since the es-band, 
which directly connects the source and drain leads, is replaced by the WS-band, where the electronic states are 
localized between two mutually inverse triangular barriers. Tunneling through es-band was already discussed in Sec. 
3. Treatment of WS-states assisted tunneling, which refers to the energy interval E < eV/2 — £, is similar to the 
preceding analysis. 

In the present case, however, large A/", explicit expressions of the Green function—matrix elements are different for 
the eigenvalue energies E n (up to exponentially small corrections given by E n w n£t3) and for E ^ E n . For the latter 
(and N = 2N + 1), we have 

vj^aiks) } ,(zl^«) exp | | < e , > , (15) 

V M {E,£)G NM {E,£) J 7rVsmh5smh<5' 

where 2 cosh 8' = eV/2 - E, <f>+ = $5 + 5v, and $ 5 / = <5'cosh<5' - sinh<5'. The use of Eq. (15) in Eq. (2) yields an 
explicit expression 




„ s 4A 2 sinh 8 sinh 8' ( 4 ^ , , 

T{E,£)k -, r-cxp --$+ , (16) 

sm 2 (irE/£ ) I [exp(5 + 6') - A 2 } 2 + A 2 (exp 5 + exp 8') > 



£ 



which provides an accurate description of the tunneling probability within WS-band except the above indicated £ 
interval and energies close to values of n£, see Fig. 6. For the eigenenergies, we have instead of Eq. (16) 

tv c ^ 4sinh(5» exp [-2 {8 n - 8jJ\ ( 4 \ 

T(n£,£) pa ^ ns exp --$„ , (17) 



sinh5Ul + ^ 2 exp(-2 ( 5„)] ^\ £ 
where $~ = - $5/ , 2cosh£„ = eV/2 + n£, and 2cosh<% = eV/2 - n£ 
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In the case of weak coupling, Eq. (17) has the same meaning for the WS-band, as Eq. (12) for the sis-band; and 
under the replacement n£ — ► E (hence, S n , S' n , <I>~ — + 5, 5',Q~ — $s — respectively) the function T(E,£) defined 
in Eq. (17) envelopes the transmission spectrum over its local maxima. The envelope is reproduced equally well by 
Eq. (9a) with the replacement of the trigonometric functions, which appear in expression (9a), by their hyperbolic 
counterparts. 

By analogy, the expression of the transmission coefficient, that follows from Eq. (16) at the energies E n = (n+l/2)£ 
(for odd Af), may be called minima envelope. Its expression is given by 



T(E,£) 



A A 2 sinh 5 sinh 5' 



[exp((5 + §') - A 2 } + A 2 (exp 5 + exp d' ) 



■ exp 



(18) 



The maxima- and minima-enveloping functions shown in Fig. 6 exhibit a striking difference in their dependence on 
energy. The reason for this can be exposed by using approximation (13) and its analogue for &s> 



3$ 5 / w (eV/2-E) 



3/2 



(19) 



the accuracy of which is illustrated in Fig. 5. The exponents of the enveloping functions are thus defined simply as 
a difference (4$~/£) and sum (4$+/£) of two FN-type exponents 4(eF/2 + Ef/ 2 /(3£ ) and 4(eF/2 - E) 3 / 2 /(3£ ). 

Notice, Eqs. (13) and (19) meet the requirements $~ "- T '" ,>Jjm " ~" " 

3/2 



and = <E>~ 



eV 



0, $+ = 2$f in = eV /(3\/2) at the mid of WS-band E = 0, 
~/3 at the top of WS-band E — eV/2. These estimates makes it easily quantifiable a 



huge quantitative difference between maxima- and minima-envelope in the mid of WS-band. 

An extremely sharp resonance structure, exhibited by the transmission spectrum of WS-states assisted tunneling, 
has the same nature as a well-known phenomenon of resonance tunneling through a barricr-well-barrier structure. 
In the given case, the barriers are of a triangular shape. The structure is totally symmetric (and the transmission 
coefficient is equal or close to unit for odd and even Af, respectively) only at the spectrum mid E = 0. The 
increase/decrease of energy strongly suppresses local maxima of electron transmission because of increasing system 
asymmetry. In contrast, because the total length of the two barriers is independent of energy, the minima envelope 
of T(E,£) depends on energy not that strongly. 

As can be concluded from Eqs. (17) and (18) and is exemplified in Fig. 6, unlike sis-assisted tunneling, the sharpness 
of resonance structure in the case of through WS-band electron transmission, to a large extent, is insensitive to the 
coupling strength. Hence, provided the symmetry conditions are met, tunneling through the mid part of a tilted 
band may serve as a nearly ideal energy filter. The anomalous sharpness and exponential decrease of the equidistant 
peaks also can be regarded as a distinctive signature of Bloch oscillations in through tilted band tunneling and related 
phenomena. 



VI. THROUGH TRAPEZOID BARRIER TUNNELING 



To wind up the discussion, we briefly consider the case of tunneling indicated in Fig. 2 by TB arrow. For energies 
outside the tilted band, E > 2 + eV/2, and for large Af, small £, the transmission probability (2) can be shown to 
take the form 



T(E,£) = 



16A 2 sinh a sinh S exp[— 2(a + 5)] 
1 + A 2 [exp(-2a) + exp(-2<5)] + A 4 exp[-2(a + 5)} 



exp 



£ 



{$6 - $a) 



(20) 



where the definition of a differs from that of 5' only by interreplacement E <-> eV/2, 2 cosh a = E — eV/2. Under 
the indicated restrictions, the energy and field dependence, given by the above relation (solid lines in Fig. 7), is in 
excellent agreement with the model exact Eq. (2), except energies close to the sis-band, where Eq. (20) predicts the 
probability of tunneling which is somewhat different from exact values of T(E,£) (not seen in the scale used). 

Equation (20) and particularly, the exponential factor exp [—4 ($5 — $ Q ) /£] looks very much differently from the 
usual WKB expression of probability to tunnel through a trapezoid barrier. Nevertheless, the WKB result can be 
retrieved from Eq. (20X- by passing to the continuous limit, in the same way, as the FN exponential factor can be 
obtained from Eq. (11) H3 Skipping quite boring calculations, we present only the final result for the WKB equivalent 
of Eq. (20) (dashed line in Fig. 7) 



T WKB {E,eV) 



16A 2 sinh 2 <5exp(-2 l 5) 
[l + ^l 2 exp(-2,5)] 2 



■ exp 



-25N ( 1 



eV 
AS 2 



(21) 
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In zero field limit, Eq. (20) transforms into 



where 2 cosh 5^ = E. Hence, without the second term in the exponent in Eq. (21), the WKB and zero-field 
expressions of T(E,£) coincide in their functional dependence on 5 and S^°\ i.e., on the imaginary electron wave 
vector within the barrier. This proves the identity of Twkb(E, eV) with Eq. (20) in zero field limit. Moreover, as 
could be expected, the WKB expression fairly well describes tunneling through a trapezoid barrier, capped by an 
(energy bounded) tilted band, up to eV values comparable with Eb w , see Fig. 7, but not the case of high voltages. 

So, Eqs. (20), (11), and (16) give an accurate explicit description of electron tunneling spacified in Fig. 2 by 
TB, FN, and WS arrows, respectively. The latter equation corresponds to a purely quantum case and therefore, it 
does not have a semi-classical analogue, as do the two former equations. For the es-band, the resonance structure of 
the transmission spectrum is shown to contain, at certain voltages, evenly spaced peaks of noncanonical WSLs the 
existence of which has been predicted in Ref. 25. 



VII. CONCLUSION 



The transmission probability, describing electron ballistic transport between two leads connected electronically via 
a single tilted band, is presented as an explicit function of electron energy, electric field parameter, thickness of the 
contact (given by a superlattice, dielectric layer, or relevant) and parameter of lead effective coupling to the contacting 
region. The derived expressions bring to light all characteristic dependencies of the tunnel event, making every point, 
one would like to know about this particular model of electric field effects on tunneling, easy for the understanding. A 
number of conclusions is made throughout the discussion and their manifold physical and experimental implications 
are illustrated in various ways. No doubt that all of them were to an extent present in numerous related studies 
but had never been exposed with the present degree of explicitncss and completeness covering all typical situations 
consistent with the model. The potential application of obtained results and methodology of their derivation is even 
more promising. Next challenge is to describe the Zener tunneling and Franz-Keldysh effect in finite, particularly, 
ultra-thin structures. Work in this direction is currently in progress. 
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Figure caption 



Fig. 1. Upper part: potential energy profile of an A/"- well supper lattice and its two lowest minibands appearing due 
to the inter-well tunneling. For a superlatice (Al^Gai-xAsAl^.Gai-x* As)tv with the well width w w = 100 A, barrier 
width u>b = 40 A, barrier hight = 0.3 eV, and electron effective mass g= 0.066m e , the lowest miniband width (for 
N >> 1) -Ebw = 3.7 meV, which is 24 times smaller, than the band gap.B 

Lower part: The lowest miniband tilted by the applied potential eV, which is smaller (up) and larger (down) than 
Ehw 

Fig. 2. Energy diagrams of lead-superlattice-lead heterostructure. Electronic spectrum of the superlattice is modeled 
by a single tilted band. At low voltages eV < i?bw (upper diagram), the electron energy can be tuned to es-band 
(resonance tunneling, RT arrow), sis-band (Fowler-Nordheim type tunneling, FN arrow), or it can be out of the 
tilted band spectrum (through trapezoid barrier tunneling, TB arrow). The high voltage case is distinct only by that 
tunneling through the mid part of tilted band is assisted by Wannier-Stark states (WS arrow). 

Fig. 3. Miniband transmission spectrum evolution under increasing voltage: from bottom to top eV = 2.00, 2.83, 
3.24, 3.46, and 3.60; N = 101. The exact level energies, i.e., solutions to Dj^(E,Af) — 0, are indicated by stars. WSLs 
with noncanonical level spacing 3£i/3, 2£i/4, §£1/5, §£1/6; and |^l/7 are labeled by open circles. 

Fig. 4. Transmission spectrum of sis-assisted tunneling through the lowest miniband of a 51-well superlattice. 

a) The miniband is tilted by the electrostatic potential difference eV of (up to down) 0.5£bw, l-5i?bw, and 2E\ >W \ 
\A\ = 0.1 (weak coupling). Open and filled circles indicate £- and 2£-spaced peaks; and those peaks, which follow 

the Airy spectrum 2 + eV/2 — E n = [3ir (n — l/4)£/2] 2/3 -£@ are indicated by squares. Dashed envelopes represent 
Eq. (12) with ££ ls replaced by E. 

b) T[E,£) is calculated from Eq. (2) (solid lines) and Eq. (11) (dashed line) for eV = 0.5i?bw and \A\ — 0.1 (weak 
coupling), 1 (intermediate coupling), and 100 (strong coupling). In both figures (a) and (b), stars indicate the values 
ofT(£f,£). 

Fig. 5. Exact (solid lines) and approximate (dashed lines) dependencies &$(E) and ^s'(E) as they are represented 
in the text. Function &s(E) is specified in Eqs. (10) (exact) and (13) (approximate); &s j (E) is represented in Eq. 
(15), and its approximate expression is given by Eq. (19). Three upper graphs correspond to sis-band. In two lower 
graphs for WS-band, the rising and lowering lines represent &g(E) and $$i(E), respectively; 3$" lm = (eV/2) 3 / 2 , 
3$" lax =(eV F ) 3/2 . Functions $ ± = $5 ± $8> are plotted by dotted lines. 

Fig. 6. Transmission spectrum of WS-states assisted tunneling through the lowest miniband of a 51-well superlattice. 
Solid and dashed oscillating lines represent exact and approximate expressions of T(E,£) given in Eqs. (2) and (16), 
respectively. Maxima envelope corresponds to Eq. (17) with n£ replaced by E, and minima envelope corresponds to 
Eq. (18). In upper part \A\ = 0.1; in the mid, |^4| = 1; and for two lower graphs, \A\ = 100. 

Fig. 7. Tunneling probability above (or below) the tilted band, calculated from Eqs. (20) and (21) for different 
voltages (as indicated), is represented by solid and dashed lines, respectively. In calculations, Af — 51, \A\ = 0.1. 
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